arXiv:1507.03807v5 [cond-mat.mes-hall] 3 Nov 2015 


Wiedemann-Franz law for magnon transport 

Kouki Nakata , 1 Pascal Simon , 2 and Daniel Loss 1 
1 Department of Physics, University of Basel, Klingelbergstrasse 82, CH-f056 Basel, Switzerland 
2 Laboratoire de Physique des Solides, CNRS UMR-8502, 

Universite Paris Sud, 91405 Orsay Cedex, France 
(Dated: November 4, 2015) 

One of the main goals of spintronics is to improve transport of information carriers and to achieve 
new functionalities with ultra-low dissipation. A most promising strategy for this holy grail is to 
use pure magnon currents created and transported in insulating magnets, in the complete absence 
of any conducting metallic elements. Here we propose a realistic solution to this fundamental chal¬ 
lenge by analyzing magnon and heat transport in insulating ferromagnetic junctions. We calculate 
all transport coefficients for magnon transport and establish Onsager relations between them. We 
theoretically discover that magnon transport injunctions has a universal behavior, i.e. is indepen¬ 
dent of material parameters, and establish a magnon analog of the celebrated Wiedemann-Franz law 
which governs charge transport at low temperatures. We calculate the Seebeck and Peltier coeffi¬ 
cients which are crucial quantities for spin caloritronics and demonstrate that they assume universal 
values in the low temperature limit. Finally, we show that our predictions are within experimental 
reach with current device and measurement technologies. 

PACS numbers: 75.30.Ds, 72.25.Mk, 85.75.-d 


I. INTRODUCTION 

The study of thermoelectric properties of materials 
started more than two centuries ago and has been receiv¬ 
ing recently enormous attention motivated by the search 
of low-power consumption technologies. One of the key 
ingredients in thermoelectrics is the Seebeck effect which 
refers to the generation of a voltage across the system 
under study when a thermal gradient is externally ap¬ 
plied to it. Conversely, the Peltier effect is tradition¬ 
ally associated with the generation of a heat flow when 
an electric current is driven across the system. Despite 
the fact that most studies focus on electric conductors, 
these fundamental concepts also extend to spin degrees 
of freedom and are at the heart of the emerging field 
of spin caloritronics which focuses on the interaction of 
spins with heat currents^— . 

A very promising strategy to look for low-dissipation 
spintronics devices is to use pure magnon currents^ which 
are collective low-energy excitations created and trans¬ 
ported in insulating magnets, therefore in the complete 
absence of any conducting metallic elements££ 

It was indeed shown that the magnetic insulator 
Y3FesOi2 (YIG) can convert a heat flow into a spin 
voltage^ which was detected by attaching a Pt thin film 
to YIG and using the inverse spin Hall effect^ to con¬ 
vert spin-wave currents into electric currents^. Such an 
observation has a natural interpretation in terms of the 
spin Seebeck effect (SSE )- 11 for magnetic insulators. In 
such materials, the SSE was shownA? — to be driven by 
magnons only which are the low-energy collective mag¬ 
netic excitation of spins. The Onsager reciprocal effect^, 
the spin Peltier effect (SPE) corresponding to the gener¬ 
ation of a magnon heat current by a spin current, has 
also been recently observed experimentally^. 

In this paper, we analyze magnon and heat currents 



FIG. 1: (Color online) Schematic representation of the fer¬ 
romagnetic insulating junction. Magnons, i.e., the quanta 
of spin-waves, carry magnetic momentum pb and heat ks- 
Those thermomagnetic transport properties are described by 
the WF law [Eq. 0j. The exchange interaction between the 
nearest neighbor spins in each FI (whose width is L ) is J and 
the magnetic field in the left (right) FI along the a-axis is de¬ 
noted by Bp r y The time reversal symmetry is broken by the 
ferromagnetic order and the magnetic field. The boundary 
spins Sr, in the left FI and Sr,, in the right FI are relevant to 
the transport of such magnons and they are weakly exchange- 
coupled with the strength J ex (<^C J). 


induced by magnetic field and temperature differences 
in insulating ferromagnetic junctions as depicted in Fig. 
□ We calculate all transport coefficients for magnon 
transport and establish Onsager relations between them. 
We predict some universal thermomagnetic relations for 
magnon transport in ferromagnetic junctions and estab- 
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lish a magnon analog of the Wiedemann-Franz (WF) 
lawil which is known to govern charge transport at low 
temperature s 18 i 19 . With respect to the SSE or SPE, we 
demonstrate that the magnon Seebeck and Peltier coef¬ 
ficients behave universally in the low temperature limit. 
We also show that these features are extremely robust 
with respect to multi-magnon effects. However, we find 
that magnon-magnon interactions can give rise to devia¬ 
tions from the Onsager reciprocity relation between the 
magnon Seebeck and Peltier coefficients. Such a mech¬ 
anism may offer an explanation for recent experiments 
by Dejene et al, 20 ' where deviations from the Onsager re¬ 
lation between the spin Seebeck and Peltier coefficients 
were observed. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the model system for a ferromagnetic junction. In 
Sec. Ill we describe the magnon transport in harmonic 
approximation driven by gradients of magnetic fields 
or heat and derive all Onsager transport coefficients 
perturbatively. There we also find the WF law and the 
universal behavior of the Seebeck and Peltier coefficients. 
In Sec. IV we discuss multi-magnon effects such as three- 
and four magnon terms. In Sec. V we give some concrete 
estimates for YIG systems. Finally, we summarize and 
give some conclusions in Sec. VI. The technical details 
are deferred to App. A-E. 

II. SYSTEM 

We consider a magnetic junction formed by two ferro¬ 
magnetic insulators (FI), as illustrated in Fig. [I] The 
temperature of the left (right) FI is T^ r ) and the cross- 
section area of the junction interface is A. Due to a 
finite overlap of the wave functions there exists in gen¬ 
eral a finite exchange interaction between the spins lo¬ 
cated at the boundaries between the two FIs. Thus, only 
the boundary spins, denoted as Sr ; and Sr r in the left 
and right FI, respectively (see Fig. |TJ), are relevant for 
magnon transport across the junction interface. The ex¬ 
change interaction between the two FIs may be described 
by the Hamiltonian2L22 Hex = -J e xJ2(r,r r ) S D ' s r r , 
where J ex > 0 is the exchange interaction, weakly cou¬ 
pling the two FIs. Assuming magnetic order along 
the magnetic field, defining the z-direction, we per¬ 
form a Holstein-Primakoff expansio n 13 ! 23 to leading or¬ 
der, S+ /r = V^Sa^ + OiS- 1 / 2 ) and Sf /r = S - a\ /r a l/r , 
where [a;, a£] = we obtain 

Hex = JexS ^ ^ ~ ar,.ka rr . k / ~t- H.c., (1) 

kx 

where k = {k x ,k y ,k z ), k' = (k' x ,k y ,k z ), k_L = (0 ,k y ,k z ), 


and the bosonic operator a £ (or r/i ) creates (annihi¬ 

lates) a boundary magnon at the right/left FI. We note 
that the Ix-momentum of magnons is not conserved at 
the sharp junction interface, whereas the perpendicular 
momentum kj^ is conserved. The tunneling Hamiltonian 
"Hex thus gives the time-evolution of the magnon num¬ 
ber operators in both FIs and generates the magnon and 
heat currents. In obtaining Eq. ©, we have assumed 
large spins S 1— and hence the 0(5'°) term in Eq. 
© indeed becomes negligible. 

Assuming cubic lattices, each of the three-dimensional 
FIs can be described by a Heisenberg spin Hamiltonian 
in the presence of a Zeeman term. The time reversal 
symmetry is broken by the assumed ferromagnetic order 
and by the magnetic field. Within the long wave-length 
approximation and in the continuum limit, the magnon 
dispersion relation in each FI reads 

oj l k {r) = 2 J Sa 2 k 2 + gn B B l(r) , (2) 

where J > 0 is the isotropic exchange interaction 
between the nearest neighbor spins in each FI, a denotes 
the lattice constant, k =| k | is the wave vector modulus, 
and Bu r ) is the magnetic field for magnons in the left 
(right) FI along the £-axis. We assume that J ex -C J 
and treat H e x perturbatively. We remark that in 
addition to the tunneling Hamiltonian given by Eq. ©, 
other 0( J ex S’)-terms actually arise from Hex due to 
the ferromagnetic order and act as effective magnetic 
fields for the boundary magnons. Including such effects, 
an effective magnetic field in Eq. 0 is introduced. 
Thus, even in the absence of an external magnetic field, 
> 0 for the boundary magnons. 

III. ONSAGER COEFFICIENTS 

The magnetic field and temperature differences defined 
by A B = B r — Bi and AT = T r — Ti, respectively, 
generate the magnon and heat current s— I m and Iq , 

where^ I m = -i( J e xS/h) ]T) k k , gpBar^kaj^ k , + H.c. 
and Iq = ~i(J ex S/h) X/k.k' w fc a ri,kar r . k / + H.c. Within 
the linear response regime, each Onsager coefficient IV 
(i,j = 1,2) and the matrix L are defined by 



A straightforward perturbative calculation in J ex ^ based 
on the Schwinger-Keldysh2L2& formalism and up to 
0(J 2 X ) gives each coefficient 
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where 


b = 


9^bB 

k B T ‘ 


(5) 


We have denoted Bi = B and T) = T for simplic¬ 
ity, and a phenomenological magnon lifetime^ r in e = 
hf3/(2r) <C 1, where /3 = (fee? 1 ) -1 , has been introduced 
and is mainly due to nonmagnetic impurity scatterings. 
Here, Li s (z) = Y^=i z n /n s is the polylogarithm function 
and Ei(— ne 2 ) = 7 + In | ne 2 \ +0(ne 2 ) the exponential 
integral, where 7 is the Euler constant. The coefficients 
are seen to satisfy the Onsager relatio n 15 ! 18 ! 24 

L 21 = T ■ L 12 . (6) 

The coefficient L 11 is identified with the magnetic 
magnon conductance G, and, in analogy to charge 
transpor t 18 ! 24 ! 26 , the thermal magnon conductance K is 
defined by K = L 22 — L 21 L 12 (L 11 ) -1 . From Eqs. (Hall - 
(I4dl) we obtain the thermomagnetic ratio K/G : charac¬ 
terizing magnon and heat transport, fts behavior is plot¬ 
ted in Fig. [5] At low temperatures^ (i.e., 1 -C fo), the 
ratio becomes linear in temperature (see also Fig. O, 


(g^/k B f [K/(GT)\ 



FIG. 2: (Color online) Plot of the ratio ( gp B /k B ) 2 [K/{GT )] 
as function of b where t — lO -10 . At low temperatures b = 
0(10), the ratio reaches the constant ‘1’ and the WF law for 
magnon transport [Eq. Q] is realized. 


K / fc R \ 2 

— = (-) T, when h/{2 r) <C k B T < gp B B. (7) 

G \g/i B J 

Therefore, in analogy to charge transport in metals^ —, 
we refer to this behavior as the Wiedemann-Franz law 
for magnon transport. The constant C analogous to the 
Lorenz number becomes 



where the role of the charge e is played by gg B . The 
magnetic Lorenz number is independent of any mate¬ 
rial parameters except the 5 -factor which is material spe¬ 
cific. Interestingly, the WF law holds in the same way 
for magnons, which are bosonic excitations, as for elec¬ 
trons which are fermions. The linear-in-T behavior can 
be traced back to the Onsager relation, Eq. and is 
independent of microscopic details^. 

Similarly, in analogy to charge transport in metal s 18 ' 24 , 
we refer to S = L 12 /L 11 as magnon Seebeck coeffi¬ 
cient (i.e., thermomagnetic power), and n = L 21 /L 11 


as magnon Peltier coefficient. The Onsager relation Eq. 
® provides the Thomson relation (i.e., Kelvin-Onsager 
relation^) 

n = TS. (9) 

At low temperatures, h/(2r) -C k B T -C gg B B, the coef¬ 
ficients reduce to 2 & 

s = n = h. (10) 

This is a remarkable result: The magnon Seebeck and 
Peltier coefficients become universal at low temperatures, 
i.e., they are completely independent of any material pa¬ 
rameters (including the 5 -factor) and are solely deter¬ 
mined by the applied magnetic field and temperature. 

Finally, we remark that at low temperatures, h/ (2r) -C 
k B T < gpL B B : each Onsager coefficient ( i,j = 1,2) 
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and the thermal magnon conductance behave in terms of 
e as 


L ij ~ lne, K ~ lne. (11) 

Thus, all coefficients show a weak logarithmic depen¬ 
dence on r, i.e., ~ Inr. In addition, Eq. © 

implies that both currents arise from terms of order 
O(Jex)- Therefore, even when an electric field is applied 
to the interface, the resulting Aharonov-Casher phased 
cannot play any significant role in the transport of 
such noncondensed magnons. Moreover, even when 
a magnetic field difference A B ^ 0 is generated, the 
noncondensed magnon current becomes essentially a dc 
one. This is in sharp contrast to the condensed magnon 
current^ which arises from the 0( J ex )-term. 

IV. MULTI-MAGNON EFFECTS 

So far we have considered the transport of essen¬ 
tially noninteracting magnons. Now we take rnulti- 
magnon effects into account. Two kinds of effects are 
considered below; the first one corresponds to a three- 
magnon splitting^ which arises due to higher order 
terms in the l/S'-expansion of the Holstein-Primakoff 
transformation^, while the second one appears because 
of magnon-magnon interaction s 21 ' 2 - £l and is due to the 
anisotropy 2 £ of the exchange interaction between neigh¬ 
boring spins in each FI. 

We begin by considering higher order terms of the 
Holstein-Primakoff transformation^, S+ = V2S\1 — 
a|a i /(4S')]a i + 0(S ~ 3 / 2 ), and thereby include the three- 
magnon splittings, a]aidi/S + H.c., into Hex- One can 
expect that each coefficient becomes smaller because such 
terms correspond to the replacement of the operator at 
with [1 — a\a,i/[AS)\ai. A straightforward calculation^ 
based again on the Schwinger-Keldysh2Zi2& formalism 
gives the simple result for the modified matrix L 2 up 
to 0(J 2 X ) 


L 2 = 


1 - 


1 fk B T^/\. , 

4^/ 2 5 V 2 JSJ Ll3/2(C 


-b\ 


) L. ( 12 ) 


Eq. m shows that, although each coefficient becomes 
smaller due to the three-magnon splittings, the modified 
matrix L 2 is still characterized by the noninteracting one 
L [Eq. (|5|)]; L 2 oc L. Consequently, the Onsager and 
the Thomson relations, Eqs. m and ©, remain satis¬ 
fied. In addition, the magnon WF law [Eq. 0] and the 
Seebeck and Peltier coefficients, [Eqs. ([Till) ], at low tem¬ 
perature remain valid. These thermomagnetic properties 
are therefore robust against the three-magnon splittings. 
One can show2£ that these properties actually hold in 
any order of the 1 /S'-expansion of the Holstein-Primakoff 
transformation. 

An anisotropi o^h^— exchange interaction between 
nearest-neighbor spins in each FI gives rise to 
magnon-magnon interactions of the typa 2 £ = 


—Jm^{ij) a i a j a i a j- The symbol (ij) indicates summa¬ 
tion over nearest-neighbor spins in each FI. The mag¬ 
nitude and the sign of the interaction J m depends on 
the anisotropy2& of the exchange interaction. We assume 
here a small anisotropy | J m | < J. A straightforward 
but tedious calculation^ gives us the modified matrix up 
to C>(J 2 x J m ) 


L 3 — 


L n +6L n L 12 
L 21 +6L 21 L 22 


(13) 


where 


Je X JmVkBTTA 


11 ( n ,,_ \ \ 2 y ; 


WVZSnWj */ 2 
x ln(i/2JS'/3aA/e), 
J 2 x J m ^/k B TrA 


(g^ B A) Li 3 / 2 (e ) 


(14a) 




SL 2L = ex “ 1V , --<?MBA 2 Li 3/2 (e- 

16V2Sir 5/2 h 2 J 5 / 2 ' ' 

x [jS'(aA ) 2 - e 2 /(2/3) + [JS'(aA ) 2 + gf i B B] 
x ln(\/2JS(3aA/e) , (14b) 


with aA = ^/5/(JS/3). These results imply the violation 
of the Onsager relation in Eq. ([6]) [and of the Thomson 
relation in Eq. @] due to magnon-magnon interactions. 
Notice that SL 11 = 0(J 2 x J m ) and SL 21 = 0(J 2 x J m ), 
while L 11 = 0(J e 2 x J°) and L 21 = 0(J 2 x J °). We re¬ 
call that the Onsager reciprocal relatio n 15 ' 35 ^ 6 could in 
principle be already violated in the noninteracting case 
since the time reversal symmetry is broken by the fer¬ 
romagnetic order and the magnetic field right from the 
outset. Still, we have microscopically found that the re¬ 
lation remains satisfied even in the presence of the three- 
magnon and higher order splitting terms of the Holstein- 
Primakoff transformation^. However, the anisotropy in¬ 
duced magnon-magnon interaction J m provides a ‘non¬ 
linearity’ 5L 21 = 0(J 2 x J m ) in terms of the perturba¬ 
tive terms (</ ex and J m ), and consequently the matrix L 3 
cannot be reduced to the form L = 0(J 2 x J^). Using 
a mean field argument or a more rigorous microscopic 
calculation^, one can show that the magnon-magnon in¬ 
teraction acts as an effective magnetic field and that the 
total magnetic field difference AB to t may be written as 
A-Btot = (1 + b m )AB with b m = 0(J m ). The term b m 
gives 6L 11 = 0(J m ) and 5L 21 = 0(J m ). The Onsager 
relation is thus violated due to the nonlinearity caused 
by the anisotropy induced magnon-magnon interaction. 
The magnitude of the effective magnetic field difference 
can be estimated by 6 m ~ 5L 21 /L 21 . 

These multi-magnon contributions, 5L 11 and SL 21 , 
generally affect also the thermomagnetic properties and 
can lead to deviations from our previous results. How¬ 
ever, at low temperatures, where the WF law, Eq. 
©, and the universality of Seebeck and Peltier coeffi¬ 
cients hold, these deviations become negligible^ because 
\S \ -C 1 for typical parameter values (see below). 

So far we have assumed bulk FIs (see Fig. [U where 
magnetic dipole-dipole interactions are negligible 3 '. 
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Such dipolar effects, however, become important in thin 
films, resulting in a modified dispersion for magnons^. 
Still, the WF law remains valid in this case too^, 
underlining the universality of this law. 

V. ESTIMATES FOR EXPERIMENTS 

The magnon currents can be experimentally measured 
by using, for instance, the method proposed in Refs. 
[a l2lll38| . Since the magnons, being moving magnetic 
dipoles with magnetic moment gg-B e zi produce electric 
fields, magnon currents can be detected by measuring 
the resulting voltage drop perpendicular to the current 
direction and magnetic field. For an estimate, we assume 
the following experimental parameter values ^ 13 ’ 29 ! 39- — : 
J = lOOmeV, J ex = lOrneV, J m = lmeV, a = lA, 
A = 3cm 2 , g = 2, r = 100ns, AT = 0.5K, and 
B = 50mT (5T) and T = 300K (0.7K) for the high 
(low^a temperature regime. Using a similar set-up as in 
Ref. [2ljj , we find that the resulting voltage drop is in the 
mV (p.V) range for high (low) temperatures. Although 
small, such values are within experimental reach and are 
actually about 10 6 ( 10 3 ) times larger than the one (^nV) 
predicted for currents in condensed magnon systems^!,. 
Alternatively, attaching a metal (e.g. Pt) to the FIs and 
using the inverse spin Hall effect^ to convert magnon 
currents into electric currents^, the magnon currents 
could also be detected&iid£ by measuring the resulting 
Hall voltage in the metal. Finally, we mention that 
the temperature difference AT can be experimentally 
produced by applying microwaves of different frequencies 
to each FI 42 ' or by local laser heatin g 39 ! 43 . 

VI. SUMMARY 

We have studied the thermomagnetic transport behav¬ 
ior of a ferromagnetic insulating junction and determined 
the Onsager coefficients in linear response regime. We 
found that at low temperatures the magnon transport 
obeys an analog of the Wiedemann-Franz law where the 
ratio of heat to magnon conductance is linear in temper¬ 
ature. Like its electronic counterpart the WF law found 
here is universal and does not depend on material param¬ 
eters except the g-factor. Quite remarkably, it exhibits 
the same linear-in-T behavior at low temperatures as the 
one for electronic transport in spite of the fact that the 
quantum-statistical properties of bosons and fermions are 
fundamentally different, in particular in the low temper¬ 
ature regime where quantum effects dominate. 

The temperature scale, however, for electrons is given 
by the Fermi temperature (~ 10 4 K for normal met¬ 
als), while it is the magnetic field for magnons (a few 
Kelvins). Obviously, these two scales are very different, 
and that is why the electronic counterpart of the WF 
law is valid typically at much higher temperatures than 
the magnonic one, but in both cases the WF law applies 
when the system temperature is low relative to its re¬ 
spective temperature scale. Moreover, we showed that 


the magnon Seebeck and Peltier coefficients become uni¬ 
versal at low temperatures. 

As an outlook we mention that it would be interesting 
to explore the regime beyond the weak junction coupling 
studied here and see if the WF law can be extended to 
such a regime as well. 

Finally, it would be interesting to test our predictions 
experimentally in candidate systems like insulating 
ferromagnets such a YIG material in the bulk or thin 
film limit. 
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Appendix A: MAGNON AND HEAT CURRENTS 

Assume cubic lattices, each three-dimensional FI is de¬ 
scribed by a Heisenberg spin Hamiltonian in the presence 
of a Zeeman term, 

RfI = — 'Y2 Si ■ J ■ S, — ^ S i; (Al) 

(ij) * 

where J denotes a diagonal 3 x 3-matrix with diag(J) = 
J{ 1 , 1 , 77 }. The exchange interaction between neighbor¬ 
ing spins in the ferromagnetic insulator is J > 0, 77 
denotes the anisotropy of the spin Hamiltonian, and 
B;( r ) = Bi( r )e~ is an applied magnetic field to the left 
(right) FI (e- denotes the unit vector along the 2 -axis). 
The symbol ( ij) indicates summation over neighboring 
spins in each FI and Sj denotes the spin of length S at 
lattice site i. Within our microscopic calculation^, we 
find that in the continuum limit, the magnon-magnon in¬ 
teraction T~l m could arise from the 77 ^ 1 anisotropic spin 
Hamiltonian Hfi as the 0(1 — 77 ) term 

Rm = — T m a 3 y dr a' l '(r)a t (r)a(r)a(r), (A2) 

where J m = —J(1 — 77 ) = 0(5°), the lattice constant 
a, and [a(r), a.t(r')] = S(r — r'). Therefore the magnon- 
magnon interaction does not influence the magnon trans¬ 
port between 77 = 1 isotropic FIs in any significant 
manne r 21 ’ 22 and we can neglect them in the isotropic 
case. Assuming the isotropic FI, within the long wave¬ 
length approximation, the magnon dispersion in each 
FI is given by Eq. © in the main text, i.e., = 

2 JSa 2 k 2 + ggB,Bi( r y 
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We then consider a magnetic junction formed by two 
ferromagnetic insulators, as illustrated in Fig. [TJ The 
temperature of the left (right) FI is T^ r ) and the cross- 
section area of the junction interface is A. Due to a 
finite overlap of the wave functions, there exists in gen¬ 
eral a finite exchange interaction between the spins lo¬ 
cated at the boundaries between the two FIs. Thus, only 
the boundary spins, denoted as Sr, and Sr r in the left 
and right FI, respectively (see Fig. [[]), are relevant for 
magnon transport across the junction interface. The ex¬ 
change interaction between the two FIs may be described 
by the Hamiltonian^^ H ex = -J e xJ2(r, ry> S D ' s r r , 
where J ex > 0 is the exchange interaction, weakly cou¬ 
pling the two FIs. Assuming magnetic order along 
the magnetic field, defining the z-direction, we per¬ 
form a Holstein-Primakoff expansio n 13 ! 23 to leading or¬ 
der, S+ /r = y/2Sa\ /r + 0(S~ 1 / 2 ) and Sf /r = S- a\ /r a l/r , 

where [a;, a£] = di tT , we obtain Eq. CD in the main text. 

The tunneling Hamiltonian "Hex gives the time- 
evolution of the magnon number operators in the FIs and 
generates the magnon and heat current s 18 ! 24 in the junc¬ 
tion. The Heisenberg equation of motion provides the 
magnon current operator 2 m and heat current operator 
Eq by 2 — 2 m = /h) ~b H- c - 

and Tq = -i(J ex S/K) ]T k X) k / w [- a r,,kar rk , + H.c. Using 
the Schwinger-Keldys h 27 ! 28 formalism and treating Hex 
perturbatively (J ex <C J), they can be evaluated up to 


o(JL) 

(E m ) = ^ Je 2 x J^ 2 a [ duj 9VB (A3a) 

k ,k f ^ 

(^Z,k,cD^r,k',cD ,cd) ’ 

wi = e / ( A3b ) 

71 k ,k’ x J 

X (^/.k^^r.k'.cj ^Z,k.oj^A,k', cj) ’ 

where L is the width of the FI (Fig. [TJ and are 

the bosonic lesser (greater) Green functions. The both 
currents arise from the 0( J 2 x )-terms; (2 m ) = 0(J 2 X ) and 
(Eq) = 0(J 2 X ). 

We phenomenologically^ introduce a life-time r of 
magnons mainly due to nonmagnetic impurity scatter¬ 
ings into Green functions (e.g., the retarded Green func¬ 
tion = [Hco — w k + ih/(2r)]~ 1 ), and regard it as a 
constant^. This gives 


Q< G A , 

^ ijkjU ;^ 7 r,k' ,u 


G> G K ,, 

^ L,k.,uj z/ r^k. ,uj 


= I 2 l^, k > I 2 

X [™("k)-n(4)], ( A4 ) 


where the Bose-distribution function n(wk) = (e^ Wk — 
1 ) _1 and /3 = l/(fcBT). Thus the currents become 


H/(2t) 


(2 m ) = — — ) V 5/ie y^[n(oJ k ) — n(oJ k ,)]-~-. 

h\ L J k [2 JSa 2 (k’ 2 - fc 2 )] 2 + [H/(2t)Y 

<x> = ft/(2r) 


(A5a) 

(A5b) 


k K [2JSa 2 (k x — k 2 )} 2 + [2/(2t)] 2 

For large 2, we can replace the sums by integrals. The summation over k! x in Eqs. (IA5al) and (IA5bl) thus becomes 

h/ ( 2 r) 


X][ n ( w k) 


„ — = —\ — » [ dx'[n(u>l) — n(ioL)] — 5 -(A 6 a) 

[2 JSa 2 (k' x 2 - fc 2 )] 2 + [ft/( 2 r )[ 2 2 tt V 2 JSa 2 J 1 K k ' V k (x' 2 - x 2 ) 2 + e* ’ 

\j 2 j So 2 J dx ^ n ^ ~ n(ul,)]n5{x' 2 - x 2 ), (A 6 b) 


L 
27r 


r 


where x 2 = 2JSa 2 f3k 2 , x' 2 = 2JSa 2 f3k x 2 , and e = 
H/3/(2t). In the last equation [Eq. (|A 6 b[) ]. we have as¬ 
sumed that r is large such that e -C 1 . 


We see that both currents [Eqs. (lA5al) and (lA5bl) ] 
are characterized by the difference of Bose-distribution 
functions n(ui). After integration over x ', Eq. (IA 6 bl) 
reduces to n(u; k ) — n(w k ). Within the linear response 


regime, this difference becomes 


n(w k ) - n(w k ) 



(e^k-l ) 2 

e^“k 

(e^“k-l ) 2 


AH, 

■AT, 


for AT = 0, 
for A B = ^ 


We have introduced the magnetic field and temperature 
differences 


AH = H r — Hi, AT = T r - T u (A 8 ) 

and have denoted Hj = H and T; = T for convenience. 
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The difference of the Bose-distribution functions has been 
expanded in powers of AT <C T and A B <C B. Eq. (IA71) 
yields 

[»(^k) - Iab=o /at = oj 

[n(u>l.) - n(w£)] |at=o /A B g^T' 

Within linear response, we can also approximate the 
magnetic field terms in all Green functions in Eq. (IA4I) 
by B. Finally, one can see that Eq. (IA9I) is responsible 
for the Onsager reciprocal relation, given in Eq. © in 
the main text. 

Thus, we see that the magnon and heat currents are 
characterized by the difference of the Bose-distribution 
functions, and the difference gives the linear responses 
[Eq. (IA7I) and (IA9I) |. and the Onsager relation, Eq. 0, 
holds accordingly. This remains in place even when three- 
magnon splittings [Eq. (IE3D ] and any higher order terms 
of the Holstein-Primakoff expansion are taken into ac¬ 
count [see Eqs. (IE1I) and (IE2I) as an example]: in any 
order of the expansion, the currents are characterized 
by the difference of the Bose-distribution functions [Eq. 
mi Consequently, the Onsager relation [Eqs. (IA9I) 
and ( 0 ] holds. 

This is changed by the anisotropy induced rnagnon- 
magnon interaction and the Onsager relation becomes 
violated, as we shall see below. 


Appendix B: ONSAGER COEFFICIENTS 

Within the linear response regime [Eq. (IA7I) ]. Eqs. 
(IA5al) and (1 A5bl) provide the Onsager coefficients I /- 7 in 
Eqs. Cal) - (l4dl) of the main text. The coefficients, L 21 and 
T 12 , are seen to satisfy the Onsager relation Eq. 0. The 
coefficient T 11 is identified with the magnetic magnon 
conductance^ G, and the thermal magnon conductance 
K is defined by K = L 22 — T 21 T 12 /T 11 . In analogy to 
charge transpor t 18 ’ 24 , this definition follows from the con¬ 
dition that the magnon current (I m ) induced by the ap¬ 
plied thermal difference AT be zero. This gives rise to an 
induced magnetic field difference ATi n( j = ATT 12 /T 41 , 
and thus to the thermal magnon current (Iq) = A'AT 
(in the absence of an applied field gradient, i.e., A B = 0 ). 

Eqs. (Hal) - (l4dl) determine the general thermomagnetic 
ratio K/G for magnon and heat transport, see Fig. [ 2 ] 
in the main text. At low temperatures, 1 -C b = 
gfj, B B / (k B T), the ratio reduces to Eq. 0 in the main 
text (see also Fig. 0. We note that the temperature 
should be still such that e <C 1 remains satisfied, i.e., 
H/{2t) -C k B T <C gg- B B. In the opposite limit, e 1, 
both currents, Eqs. (IA5al) and (IA5bl) . are seen to be 
exponentially vanishing as T —> 0. 

The Onsager relation [Eq. ©] provides the Thomson 
relation (i.e., the Kelvin-Onsager relation^) in Eq. 0 of 
the main text. At low temperatures, H/(2t) <C k B T <C 
gg B B, the magnon Seebeck and Peltier coefficients re¬ 


duce to Eq. (fTUl) in the main text (see Fig. 0, and, quite 
strikingly, become universal. 

We remark that at low temperatures, each Onsager 
coefficient T *- 7 and the thermal magnon conductance be¬ 
have in terms of e as in Eq. (HID in the main text. Thus, 
all coefficients show a weak logarithmic dependence on 
r, i.e., T 1 - 7 , K ~ lnr. 


gju B S/(k B b) 



FIG. 3: (Color online) Plot of magnon Seebeck coefficient 
S as function of b where e = 10 -10 . At low temperatures 
b = O( 10 2 ), the rescaled coefficient gpB*S/(fcB&) approaches 
unity asymptotically and realizes the universal relation [Eq. 

G0[. 


Appendix C: YIG THIN FILM 

So far we have assumed that both FIs are of bulk shape 
such that surface effects due to magnetic dipole-dipole 
interactions are negligible 37 . However, it is interesting 
to consider also a thin film geometry since such systems 
are of great experimental interest, in particular for YIG 
films^. 

We show now that the WF law also holds for a thin 
film, where now the length of the junction in x-direction 
L (see Fig. 0 is short compared to all other dimen¬ 
sions. Due to the dipole-dipole interaction, the magnon 
dispersion relation in each FI changes and becomes^ 
uj 1 ^ = D(k 2 — k^) 2 + g^BB^ r ), where fc m ~ 10 4 /cm 
for e.g. YIG thin films^. The parameter D is due 
to the long-range dipole-dipole interaction as well as 
the exchange interaction between the nearest-neighbor 
spins. The main contribution comes from k ~ fc m where 
« ADk‘$ a {k—k m ) 2 +gg,BB l ( r ) for k ~ fc m . As we have 
seen in the main text, the WF law is realized in strong 
magnetic fields T;( r ) ~ a few T (leading to a large en¬ 
ergy gap). Repeating the same perturbative calculation 
as before, but under the restriction that only the low¬ 
est magnon subband in ^-direction (k x = k' x = 0 ) needs 
















to be taken into account due to finite-size quantization coefficients in thin films become 
(which actually simplifies the calculation), the Onsager 

I 


L 11 
L 12 
L 21 
A 22 


tA I 

( Jen.^Q' 

47rh 2 Dkl ' 

V L 

tA | 


4tt WDk^ ' 

V L 

tA | 


4tt hADk^ ’ 

V L 

tA I 

( JeyiS(l 

AirhADk^ ’ 

V L 


( 5 /x B ) Li 0 (e b ) 


\ 2 r 

J gg B&B Lii(e _b ) + 6Li 0 (e _b ) 

) g^k B T Lii(e _b ) + 6Li 0 (e _b ) 

) : 


klT 


2Li 2 (e _f> ) + 2fcLii(e - °) + ^Lio^" 6 ) 


(Cla) 

(Clb) 

(Clc) 

(Cld) 


The Onsager relation, L 21 = T ■ L 12 , again holds. 

In particular, the ratio K/G reduces to^ 

K ^ / k B \ 2 

— = (-) T, when k B T < gg, B B. (C2) 

G \gg, B J 

Again, the ratio is linear in temperature, and, quite re¬ 
markably, again with a universal prefactor (in particular 
independent of the dipole-dipole interaction D). Thus, 
we see that the WF law holds also for thin films and in 
the presence of dipole-dipole interactions. Thus, we con¬ 
clude that the linear-in-T behavior is extremely robust 
against microscopic details. 

We note that each Onsager coefficient itself [Eqs. 
(ICla|l - (|Cld|l | drastically changes in terms of temperature 
dependence as well as magnon lifetime etc. as compared 
to the bulk results (see main text). However, these dif¬ 
ferences all cancel out in the ratio K/G and the linear- 
in-temperature behavior remains. 

Finally, at low temperatures k B T -C gg B B, the 
magnon Seebeck and Peltier coefficients for thin films 
are reduced to S = B/T and n = B as for the bulk case 
given in the main text. 


up to 0(J 2 x J m ) 


,-T x (JexS ) 2 j> 2 fk B Ty 

{m> 2 ^/ 2£2 \2JSJ 


Li3/2(e b ) 


k ,/c' 


X R e(GI,k, u Gl,k,w3r,k '+ Gl,k,uGl,k.u>G//k',u 

i r>< n r /->a, , x->< ^-»a \ 

' yl,let,LU^rM' ,UJ ' ,UJ ) 

- N r), 


duigfj, B 


(Dl) 


where c; a ( r /</>) is the bosonic advanced (re¬ 
tarded/lesser/greater) Green function. Thus the 
magnon current (X m ) = 0(Jn/"^) arises as the nonlinear¬ 
ity in terms of the perturbative terms (J^ and J ex ). 
After some straightforward manipulations [Eq. (IA6bl) [. 
Eq. (1D1I) is reduced to the form 


(I m ) = + J/ a N r ) 


(D2) 


+ dk C^Af 1 - J r m Af r )[n( 4) - «K)] 


A B, 


where Ci( 2 ) is a AH- and AT-independent constant, and 
A fd r ) defined by 


Appendix D: MAGNON-MAGNON 
INTERACTIONS 

Assuming an anisotropio 21 ! 22 ' 34 exchange interaction 
among the nearest neighboring spins in each FI [Eq. 
mi magnon-magnon interactions in Eq. (IA2I) may 
arise from such an anisotropic Heisenberg spin Hamil¬ 
tonian. The sign and the magnitude of the magnon- 

l(r) 

magnon interaction Jm ' depends on the anisotropy in 
the left (right) FI and it is assumed to be small | J^ \ 
•C J to be treated perturbatively; in Eqs. (I14al) and 
(I14bll . they are set J l m = = J m for simplicity. 

Using the Schwinger-Keldysh formalis m 27 i 28 , we find 


Ul = 4^(^l) U3 / 2(e b) ’ (D3a) 

M r = M l + O(AB) + O(AT), (D3b) 

is obtained by integrating the Bose-distribution func¬ 
tion n(w^ r ’ > ) over the wavevector. Eq. (1D2I) indicates 
that the perturbation of the magnon-magnon interac¬ 
tion Jni r) work as the magnetic field difference A B ex¬ 
pansion. A mean field analysis actually shows that the 
magnon-magnon interaction works as an effective mag¬ 
netic field. The Ci-term arises from Q\ k 0J S[ k Jd\. k , u in 
Eq. (ID1I) and the C 2 -term from the rest. Since the dif¬ 
ference of the Bose-distribution functions n(w k ) — n(w k ) 
indeed gives the responses to A B and AT by itself [Eq. 
mi the C 2 -term becomes irrelevant in the linear re¬ 
sponse regime; the response to the temperature difference 































9 


AT actually arises also in the magnon-magnon interac¬ 
tions (i.e., from the C 2 -term), but it is accompanied with 
the magnetic field difference as such A B ■ AT. Therefore 
any linear responses to the temperature difference in the 
magnon-magnon interactions cannot be found in the cor¬ 
responding magnon current [Eq. or the modified 

matrix L 3 [Eq. (TT51) |. We remark that the same remarks 
apply to the difference of magnon-magnon interactions 
A J m = — J l m . The response to A J m actually arises 

from the C 2 -term [Eqs. (ID2[) . (ID3b[) . and (|A7[) [. but it is 
accompanied with nonlinear responses such as A B ■ AT 
and (AS) 2 . Therefore, any responses to AJ m cannot 
affect the modified matrix I 3 [Eq. (1 1j [) j. 

Finally, focusing on the linear response regime, the 
magnon current due to the magnon-magnon interactions 
becomes 


<I m ) = + J r m N r )AB. (D4) 

This gives SL 11 in Eq. (I14al) . where they are set ./( n = 
J ^ = J m . Using SL 11 , the constant is given by C\ = 

-6L^/( JIN 1 + rjf). 

The main mechanism behind the generation of the co¬ 
efficient SL 21 due to the magnon-magnon interactions re¬ 
mains in place. The heat current ( Iq ) = can be 

evaluated in the same way and it is indeed expressed in 
the same form as Eq. m in terms of the Green func¬ 
tions up to 0 (J 2 X J m ) 


(Iq) 


(■ J ex S) 2 J l m a 2 
2n 5 / 2 L 2 


' k B T\3/2 

Jus) 


Li3/2(e 


—b 


)£ 

KK 


dujuj l k 


X R e(Gi'k >u Gl t k tU1 Gr t k> lU + ^/ r ,k,w^ r ,k,tj^k',w 

1 n< n't nSb , n< n a n a \ 

+ ^kjW^k^^k'.w ^Z,k,o;^Z,k,cj^r,k',uJ 

- (Zor). (D5) 


After some straightforward manipulations with the help 
of Eq. (IA 6 bl) , it finally reduces to 


( 1 Q ) = \c[{JlN l + r m N r ) 


(D 6 ) 


dk C' 2 (J l m M l - j;Ar)K4) - n(£4)]l A B, 


where Cq 2 ) is the corresponding AB- and AT- 
independent constant. This result has the same quali¬ 
tative form as the magnon current Eq. (ID2I) . Within 
linear response, it becomes [Eq. (1A7I) [ 


difference AB to t may be written as AB to t = (l + & m )Af?, 
where b m = 0(J l m N l + AT r ) is the contribution of 
such magnon-magnon interactions. This contribution b m 
gives rise to the C^-term [Eqs. (ID4I) and (ID7I) . and 
Eq. (ID3al) [. and leads to SL 11 = 0(J l m N l + J^ n Af r ) and 
<5L 21 = + J^ n M r ). Thus, the Onsager relation is 

violated by the magnon-magnon interactions. The mag¬ 
nitude of the effective magnetic field difference can be 
estimated by b m ~ SL 21 /L 21 . 

The contributions arising from magnon-magnon inter¬ 
actions, SL 11 and SL 21 , generally affect also the ther- 
momagnetic properties of magnon and heat transport. 
However, at low temperatures, H/(2t) -C k B T <C g^i B B 7 
the ratios reduce to SL 11 /L 11 ~ (3J m r /2H)(k B T / S J ) 3 / 2 
and similarly for SL 21 /L 21 . For typical parameter values 
(see main text) we find that e.g. SL^/L 11 ~ 0.1. 
Therefore, L 3 = L at low temperatures. Thus, the 
Onsager and the Thomson relations, Eqs. © and ©, 
the WF law for magnon transport, Eq. ©, and the ther- 
momagnetic properties of magnon Seebeck and Peltier 
coefficients, Eq. CUD. still hold at low temperatures. 


Appendix E: THREE-MAGNON SPLITTINGS 


Until now, we have essentially used the linearized^ 
Holstein-Primakoff transformation Sf = \/2S[\ — 

a\ai/ (25')] 1 / 2 ai = \/25a, + 0(S)~ 1 / 2 . Now, we consider 
the higher term 0{S)~ 1 / 2 , S) + = \/25[l — ajai/(4S)]aj+ 
0(S)~ 3/2 , and include the three-magnon splittings, 
a\aiai/S + H.c., into the Hamiltonian "H ex to be treated 
perturbatively. We remark that when S —> 00 (i.e. in 
the classical limit), the three-magnon splittings cease to 
work. In that sense, the three-magnon splitting can be 
regarded as a quantum effect^. 

Following the same procedure as before with the ap¬ 
proximation Eq. C£zD , the magnon current in linear re¬ 
sponse regime becomes in leading order 



which is of order 0(J 2 X ). The heat current can be eval¬ 
uated in the same way and becomes up to 0( J 2 X ) 


(I Q ) = C[(J l m AT l + J r m Af r )AB. (D7) 

This gives SL 21 in Eq. (I14bl) . where J l m = J^= J m . The 
Onsager relation is violated by SL 21 due to the magnon- 
magnon interactions. Using SL 21 , the constant is given 
by C[ = —SL 21 /(J l m N l + J r m AU). 

These results can be understood as follows; the 
magnon-magnon interaction Jm ’ 1 acts as an effective 
magnetic field, and consequently the total magnetic field 


(Iq) = 


1 - 


1 /k B T\3/2 

47T 3 / 2 S\2JS) 3/2 (C 


-b\ 


x (!’■ £”) ( a A t B ) . 


Thus, the modified matrix L 2 is given by [Eq. ©] 


L, = 


1 - 


1 fk B T^/2 T f __ b< 

47t 3 / 2 5 V 2JS / 13/2 ^ ^ 


L. 


(E2) 


(E3) 
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This means 

L 2 oc L (E4a) 

= L , when T —> 0. (E4b) 


© and ©] hold. In addition, the WF law for magnon 
transport [Eq. (0)] and the thermomagnetic properties 
[Eq. (1101) ] at low temperatures remain valid even in the 
presence of three-magnon splittings. 


Therefore, the Onsager and the Thomson relations [Eqs. 
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